function [a,da,dda] = func_InteVV(dt,a,da,dda,Q,M,K,C,MC)
%myFun - Description
%
% Syntax: [a,da,dda] = func_InteVV(a,da,dda)
%
% Long description
    
 % 核心递推公式
    dahl=da+0.5*dda*dt;
    al=a+dahl*dt;
    
    % 在t+1时刻上求解加速度
    QQ=Q-K*al-C*dahl; % 如果是时间变化的外力，此处需修改
    ddal=MC\QQ;
    
    % 在t+1时刻上求速度
    dal=dahl+0.5*ddal*dt;
    
    % 时间推进
    a=al;
    da=dal;
    dda=ddal;

end